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On Exclusive and Inclusive Rare B Decays: 
CKM Phenomenology and New Physics Reach 
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We report recent results in rare B decays. Our focus will be on 6 ^ d"f and b st^i~ transitions. We discuss their 
impact on the CKM phenomenology and their role within our search for new physics. In particular, we analyse the impact 
of a recent lattice QCD estimate of the B — > K* form factor at zero recoil. We also briefly discuss the presently available 
optimal theoretical tools for the inclusive and exclusive modes. 



1 Introduction 

Rare B decays, as flavour changing neutral cur- 
rent processes (FGNG), are the most important tools 
within our (indirect) search for new physics in the 
present era of B factories (for a recent review see [1]). 
Among them, the inclusive b s"f mode is still the 
most prominent, because it has already measured by 
several independent experiments [2-6]. The present 
world average is [7] : 



BiB Xsj) = (3.34 ± 0.38) x 10" 



(1) 



The present next-to-leading -logarithmic (NLL) QGD 
prediction, based on the original QGD calculations of 
several groups ([8-11]), has an additional charm mass 
renormalization scheme ambiguity as was first noticed 
in [12] . The AI S scheme is used in the most recent up- 
dates [12,13]. This choice is guided by the experience 
gained from many higher order calculations where the 
mass is dominantly off-shell and leads to the following 
theoretical prediction for the 'total' branching ratio: 



B{B Xsj) = (3.70 ± 0.30) x 10" 



(2) 



However, the renormalization scheme for tUc is an 
NNLL issue, and one should regard the preference 
for the MS scheme in comparison with the pole mass 
scheme (which was used within all previous NLL pre- 
dictions) just as an educated guess of the presently 
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unknown NNLL corrections. Therefore, one could also 
argue for a slightly larger theoretical error in (2). A 
complete resolution of this problem can be achieved by 
a NNLL calculation which is presently under study. 

The stringent bounds obtained from the B X^j 
mode on various non-standard scenarios (see e.g. [14- 
18]) are a clear example of the importance of clean 
FGNG observables in discriminating new-physics mod- 
els. The exclusive modes, however, often have large 
uncertainties due to the hadronic form factors and are 
not clean enough to disentangle possible new physics 
effects from hadronic uncertainties, but they can serve 
as important QGD tests. Exceptions are ratios of ex- 
clusive quantities like asymmetries in which a large 
part of the hadronic uncertainties cancel out. 

Regarding the corresponding exclusive mode B 
K*^, quite recently a preliminary lattice determina- 
tion of the B — > K* form factor at zero recoil [19] 
was presented which is in perfect agreement with the 
previous indirect determinations. We analyse the im- 
pact of this new result on the GKM phenomenology 
in section 3.2. 

Besides the 6 — > S7 mode, also the b sl^t" tran- 
sitions are already accessible at B factories [20-22], 
inclusively and exclusively. The 6 — > s£'^£~ mode rep- 
resents new sources of interesting observables, particu- 
larly kinematic observables such as the invariant dilep- 
ton mass spectrum and the forward-backward (FB) 
asymmetry. Rare B decays are also relevant to the 
GKM phenomenology; the 6 d7 is especially impor- 
tant in this respect. In the following, we will focus on 
the latter two rare modes and also briefly discuss the 
theoretical tools available for the analysis of exclusive 
and inclusive channels. 
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2 Theoretical Tools 

The effective field theory approach serves as a the- 
oretical framework for both inclusive and exclusive 
modes. The standard method of the operator product 
expansion (OPE) allows for a separation of the B me- 
son decay amplitude into two distinct parts, the long- 
distance contributions contained in the operator ma- 
trix elements and the short-distance physics described 
by the so-called Wilson coefficients. The W boson and 
the top quark with mass larger than the factorization 
scale are integrated out, i.e. removed from the the- 
ory as dynamical fields. The effective Hamiltonian for 
radiative and semileptonic 6 — > s/rf transitions in the 
SM can be written as 



K 



eff 



4Gf 

"71" 



10 



A*^ao, + A«^a((!)i-on 



i=l 



(3) 



where Oi{^.) are dimension-six operators at the scale 

fi ^ O(TOfc); Ci{i-i) arc the corresponding Wilson co- 
efficients. Clearly, only in the sum of Wilson coeffi- 
cients and operators, within the observable H, does 
the scale dependence cancels out. Gp denotes the 
Fermi coupling constant and the explicit CKM fac- 



tors are A* = 
relations A^ 



V,,V,* and A^ 



V^^V*^. The unitarity 



A" were already used by us in 



(3). 

The operators can be chosen as (we only write the 
most relevant ones): 



02 


= {SLlnCL){cLl^hL) , 


(4) 




= {SLlij,UL){uLj''bL) , 


(5) 






(6) 




= l/gsm^{sLCT^'"T''hR)Gl, , 


(7) 


Os 


= e^/9l{sLl^hL)Y.{hn), 


(8) 






(9) 



where the subscripts L and R refer to left- and right- 
handed components of the fermion fields. In 6 ^ s 

transitions the contributions proportional to A" are 
rather small, while in 6 ^ d decays A^ is of the same 
order as A^ and they play an important role in CP 
and isospin asymmetries. The operators Og and Oiq 
only occur in the semileptonic b s/dl'^fr modes. 

While the Wilson coefficients Ci(/x) enter both in- 
clusive and exclusive processes and can be calcu- 
lated with perturbative methods, the calculational ap- 
proaches to the matrix elements of the operators differ 
in both cases. Within inclusive modes, one can use 
quark-hadron duality in order to derive a well-defined 



heavy mass expansion of the decay rates in powers of 
^Qco/'nib (HME). In particular, it turns out that the 
decay width of the B is well approximated by 

the partonic decay rate, which can be calculated in 
renormalization-group-improved perturbation theory: 

T{B X,7) = r(6 ^ X|'"''*°"7) + A"°"P'=''*- (10) 

Non-perturbative effects, A"°"^'^'"*- , are suppressed by 
inverse powers of ruh and are well under control thanks 
to the Heavy Mass Expansion (HME); they can be fur- 
ther estimated through the application of the Heavy 
Quark Effective Theory (HQET). In exclusive pro- 
cesses, however, one cannot rely on quark-hadron du- 
ality and has to face the difficult task of estimating ma- 
trix elements between meson states. A promising ap- 
proach is the method of QCD-improved factorization 
that has recently been systemized for non-leptonic de- 
cays in the heavy quark limit. This method allows for 
a perturbative calculation of QCD corrections to naive 
factorization and is the basis for the up-to-date predic- 
tions for exclusive rare B decays. However, within this 
approach, a general, quantitative method to estimate 
the important I /mi, corrections to the heavy quark 
limit is missing. More recently, a more general quan- 
tum field theoretical framework for the QCD-improved 
factorization was proposed - known under the name of 
Soft ColUnear Effective Theory (SCET). 

2.1 Inclusive Modes 

In contrast to the exclusive rare B decays, the inclu- 
sive ones are theoretically clean observables and dom- 
inated by the partonic contributions. Bound-state ef- 
fects of the final states are eliminated by averaging 
over a specific sum of hadronic states. Moreover, also 
long-distance effects of the initial state are accounted 
for, through the heavy mass expansion in which the 
inclusive decay rate of a heavy B meson is calculated 
using an expansion in inverse powers of the b quark 
mass. 

The optical theorem relates the inclusive decay rate of 
a hadron Hf, to the imaginary part of certain forward 
scattering amplitudes 



T{Hb ^ X) 



1 



2mH, 



■Im {Hb\T\ Hb) . 



(11) 



where the transition operator T is given by T = 
i J d'^xT[Heff{x)Heff{0)]. It is then possible to con- 
struct an OPE of the operator T, which gets expressed 
as a series of local operators - suppressed by powers of 
the b quark mass and written in terms of the b quark 
field: 



OPE 



1 

nib 



Oo 



1 

nib 



1 



r02 + ...) 



(12) 
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This construction is based on the parton-hadron du- 
ahty, using the facts that the sum is done over all 
exclusive final states and that the energy release in 
the decay is large with respect to the QCD scale, 
Aqcd < rnh. With the help of the HQET, namely 
the new heavy-quark spin-flavour symmetries arising 
in the heavy quark limit rUb ~^ oo, the hadronic ma- 
trix elements within the OPE, {Hi, \ Oi \ Hb), can 
be further simplified. The crucial observations within 
this well-defined procedure are the following: the free 
quark model turns out to be the first term in the con- 
structed expansion in powers of I /nib and therefore 
the dominant contribution. This contribution can be 
calculated in perturbative QCD. Second, in the appli- 
cations to inclusive rare B decays one finds no correc- 
tion of order l/rUb to the free quark model approxi- 
mation, and the corrections to the partonic decay rate 
start with 1 /ml only. The latter fact implies a rather 
small numerical impact of the non-perturbative cor- 
rections to the decay rate of inclusive modes. 

The l/ml corrections correspond to the OPE for 
T(0j07). There are additional non-perturbative ef- 
fects if one also takes into account the operator 02 • 
They can be analysed in a model-independent way and 
scale with l/m^. Due to small coefficients in the ex- 
pansion also their impact is very small [23-25]. 

2.2 Exclusive Modes 

The naive approach to the computation of exclusive 
amplitudes consists in writing the amplitude A ~ 
Ci{^b){Oi{nb)) and parametrizing {Oi{nb)) in terms of 
form factors. A substantial improvement is obtained 
through the QCD-improved factorization [26,27] and 
SCET [28-33] approaches. 

Let us consider processes involving the decay of a 
heavy meson into fast moving light particles {B 
^ev, B {p,K*)'-f, B Kn, ...) and indicate with 
Q ~ 0{mb) their typical large energy. The idea is to 
isolate all the relevant degrees of freedom necessary to 
correctly describe the infrared structure of QCD below 
the scale Q and associate independent fields to each 
of them. It is possible to identify two distinct pertur- 
bative modes, called hard (j>^ ~ Q^) and semi-hard 
{p^ ~ AqcdQ)- These modes are produced, for in- 
stance, in interactions of energetic light particles with 
the heavy quark and the i?-mcson spectator, respec- 
tively. These two modes do not appear in the initial 
and final states and, therefore, have to be integrated 
out. We do not wish to entertain here a comprehen- 
sive discussion of the technicalities involved in this 
step. It will suffice to say that the resulting theory 
(also called SCET// in the literature) contains only 
non-perturbative degrees of freedom with virtualities 
^(Aqcd) ^■^'i ^^^^ hard and semi-hard modes are re- 



flected in the coefficient functions in front of the op- 
erators of that (SCET//) theory. We note that these 
coefficients depend, in general, on energies of order Q 
and Aqcd- Moreover, the hierarchy Aqcd ^ Q allows 
for an expansion in the small parameter A = Aqcd/Q- 

Given a process, one has to construct the most gen- 
eral set of (SCET//) operators at a given order in A, 
and show that all the possible gluon exchanges can 
be reabsorbed, at all orders in perturbation theory, 
into form factors and meson light-cone wave func- 
tions. The resulting amplitude is a convolution of 
these non-perturbative universal objects with the co- 
efficient functions encoding the contribution of hard 
and semi-hard modes. Questions regarding the conver- 
gence of these convolution integrals can be addressed 
using symmetries, power counting and dimensional 
analysis (a discussion on this point is presented in 
Ref. [33]). 

The few form factors that describe the transition B 
M (where M denotes a pseudo-scalar or vector meson) 
can be written as [34]: 

Ff^^ =Ci^''^^ + cl,B^Ti^cPM+o(—) (13) 

where ^^^^ is the so-called non-factorizable (or soft) 
contribution to the form factors (actually there is one 
soft form factor for the decay into pseudoscalar meson 
and two for the decay into vector mesons); <j)BM are 
the B and M meson light-cone wave functions; Ci are 
Wilson coefficients that depend on hard scales; and Tj 
are perturbative hard scattering kernels generated by 
integrating out hard and semi-hard modes. In Ref. [35] 
the factorization formula Eq. (13) has been proved at 
all orders in perturbation theory and at leading order 
in AQco/mb, using SCET techniques.^ The strength 
of Eq. (13) is that it allows us to express several inde- 
pendent QCD form factors in terms of only one soft 
form factor (two in the case of vector mesons) and 
moments of the light-cone wave functions of the light 
pseudo-scalar (vector) and B mesons. 

Let us now briefly discuss the form of factorization for 
the decays B ^ V-f (with V = K*, p). At leading 
order, only the operator contributes and its matrix 
element between meson states is given by an expres- 
sion similar to (13). The choice of using either the full 
QCD form factor T^^^ or the soft one clearly is a 
matter of taste (note that non-perturbative methods, 
such as lattice-QCD and light-cone QCD sum rules, 
give only informations on the full QCD form factors 
and not on the soft contributions alone). The advan- 
tage of the QCD-improved factorization approach is 

^We note that a discussion of the convergence of the convolution 
integrals is still missing. 
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evident in the computation of the next-to-leading or- 
der (in tts) corrections. In fact, one can show that the 
matrix elements of the operators O2 and Og, which are 
expected to contribute at this order, are given by the 
matrix element of O7 times a computable hard scat- 
tering kernel. Moreover, spectator interactions can 
be computed and are given by convolutions involving 
the light-cone wave functions of the B and V mesons. 
It must be mentioned that light-cone wave functions 
of pseudo-scalar and vector mesons have been deeply 
studied using light-cone QCD sum rules methods [36- 
39]. On the other hand, not much is known about 
the B meson light-cone distribution amplitude, whose 
first negative moment enters the factorized amplitude 
at NLO. Since this moment enters the factorized ex- 
pression for the _B ^ 7 form factor as well, it might 
be possible to extract its value from measurements of 
decays like B jeu, if it can be shown that power 
corrections are under control [40]. 

Finally, let us stress that a breakdown of factorization 
is expected at order K,^ao/mb [27,41,42]. In Ref. [41], 
in particular, the authors have shown that in the anal- 
ysis of B — > K*^ decays at subleading order an in- 
frared divergence is encoTintered in the matrix element 
of Og. Nevertheless, some very specific power correc- 
tions might still be computable. Indeed, this is the 
case for the annihilation and weak exchange ampli- 
tudes in S ^ p7 at the one-loop level. 



3 6 — > (^7 Transitions 



3.1 B ^ Xd7 

Most of the theoretical improvements on the per- 
turbative contributions and the power corrections in 
l/m\ and carried out in the context of the 

decay B Xg'y, can straightforwardly be adapted 
to the decay B X^^y] thus, the NLL-improved 
decay rate for B — > ^^7 decay has much reduced 
theoretical uncertainty [43]. But as = V^^y*^ for 

^tb^td and 



6 ^ 6^7 is not small with respect to 
A^ = V^f^V*^, one also has to take into account the 
operators proportional to AJ| and, moreover, the long- 
distance contributions from the intermediate u-quark 
in the penguin loops might be important. However, 
there are three soft arguments that indicate a small 
impact of these non-perturbative contributions: first, 
one can derive a model-independent suppression fac- 
tor Aqcd/™?) within these long-distance contributions 
[44]. Second, model calcidations, based on vector 
meson dominance, also suggest this conclusion [45]. 
Furthermore, estimates of the long-distance contri- 
butions in exclusive decays B ^ and B luj 
in the light-cone sum rule approach do not exceed 
15% [46]. Finally, it must be stressed that there 



is no spurious enhancement of the form log(m„/^h) 
in the perturbative contribution, as was shown in 
[9,47]. All these observations exclude very large long- 
distance intermediate w-quark contributions in the 
decay B X^j- Nevertheless, the theoretical status 
of the decay B X^^y is not as clean as that of 
B ^ X,7. 

While the b ^ s transitions like B Xsj have no 
relevant impact on the CKM phenomenology because 
of the flatness of the corresponding unitarity trian- 
gle (for example: Vts cannot be further constrained 
by the B — > ^^7 measurement if the unitarity con- 
straint is not used in the theoretical prediction, be- 
cause the charm contribution is twice as large as the 
top contribution), b ^ d transitions give important 
complcinicintary information on the unitarity triangle, 
which is also tested by the measurements of Vub/Vcb, 
AMb^, and AMs^/AAfs^. Thus, a future measure- 
ment of the B Xd"f decay rate will help to reduce 
the currently allowed region of the CKM Wolfenstein 
parameters p and rj significantly. 

Regarding new physics, the branching ratio of B — > 
Xdj might be of interest, because its CKM suppres- 
sion by the factor ]Vtd|^/|14s|^ in the SM may not be 
true in extended models. We also emphasize that in 
the ratio 



B{B ^ Xrf7) 
B{B^X,^y 



(14) 



a good part of the theoretical uncertainties cancel out. 

It is therefore of particular interest for CKM phe- 
nomenology and for the new physics search. 

A measurement of the B X^^y is rather difficult but 
perhaps within the reach of the high-luminosity B fac- 
tories. Such a measurement will rely on high statistics 
and on powerful methods for the kaon-pion discrimi- 
nation. At present only upper bounds on correspond- 
ing exclusive modes are available (see next section). 

The direct normalized CP asymmetry of the inclusive 
decay modes represent another interesting observable: 



acp 



T{B X,/,7) - T{B X,/,-7) 
T{B^X,,d^) + T{B^X„^y 



(15) 



CLEO has already presented a measurement of the 
CP asymmetry in the inclusive decay B Xs^, ac- 
tually a measurement of a weighted sum, acp = 
QMbacp{B X,7) + 0.02acp(S ^ Xd^i) [48], 
which already excludes very large effects. The same 
conclusion can be deduced from the measurements of 
the CP asymmetry in the exclusive modes (see next 
section). 
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Theoretical NLL QCD predictions of the normalized 
CP asymmetries of the inclusive channels ([43,49]) 
within the SM can be expressed by the approximate 
formulae (see [50]): 



acp{B 
acp{B 



Xal) 



0.334 X Im [e^] 
0.334 X Im [ej 



+0.6% . 
-16%, 



(16) 



where 



-x'^{p-i'n), cd = 



v*v, 



dVub 



1 



- ITj 



Numerically, the best-fit values of the CKM parame- 
ters are used. The two CP asymmetries are connected 
by the relative factor X'^ {{\ — pY + rf) . Moreover, the 
small SM prediction for the CP asymmetry in the de- 
cay B Xg') is a result of three suppression factors. 
There is an as factor needed in order to have a strong 
phase; moreover, there is a CKM suppression of order 
and there is a GIM suppression of order (mc/m(,)^, 
reflecting the fact that in the limit rric = m„ any CP 
asymmetry in the SM would vanish. 

It will be rather difficult to make an inclusive mea- 
surement of the CP asymmetry in the b ^ d channel. 
However, based on CKM unitarity, one can derive the 
following U-spin relation between the un-normalized 
CP asymmetries [51]: 



AT{B ^ X,7) + AT{B ^ X^j) = 0. 



(17) 



Within the inclusive channels, one can rely on parton- 
hadron duality and can actually compute the U-spin 
breaking by keeping a non-vanishing strange quark 
mass [52]. Going beyond the leading partonic con- 
tribution, one can further check if the large suppres- 
sion factor from the U-spin breaking is still effective, 
in addition to the natural suppression factors already 
present in the corresponding branching ratios [53] ; this 
finally leads to the SM prediction 



AT{B ^ X,j) + AT{B ^ Xdj) = 1 x 10" 



(18) 



This prediction provides a very clean SM test, whether 
generic new CP phases are active or not. Any signif- 
icant deviation from the estimate would be a direct 
hint of non-CKM contributions to CP violation. 



3.2 B 



PI 



In the analysis of exclusive B — > decays (with V = 
K*, p, uj) we will construct the various observables 
in terms of the CP- averaged quantities - which are 
much easier to measure than the individual channels 
- unless otherwise stated. In the NLL approximation. 



this procedure is equivalent to defining two distinct 
observables for the charge-conjugate modes and then 
perform the average. The ratios R{pj/ K*j) are given 

by 



R\pi/K*^) = - 



td 



{Ml 



Vts 

Vtd 



en 



{M%-Ml,f 



Ml,f 



Ai?±), 
\- AR°) , 



where C = C'[(0)/Cf * (0), and ^1(0) and Cf*(0)) are 
the form factors at = in the effective heavy 
quark theory for the decays B — > p{K*)^ [54]. There 
are several estimates of the qiiantity C, in the present 
literature coming from light-cone QCD sum rules 
(LCSR) [55], hybrid LCSR [56], improved LCSR [57] 
and quark models [58]. In the numerical analysis we 
adopt the value C = 0.76 ± 0.10; the central value is 
taken from the LCSR approach while the error is in- 
creased in order to accommodate all the other determi- 
nations. The quantities (1 -|- Ai?^''') entail the explicit 
0{as) corrections as well as the power-suppressed an- 
nihilation contributions proportional to A^. The lat- 
ter contribution, in particular, is effective only for 
charged B decays (weak annihilation); in fact, W- 
exchange amplitudes are smaller because of the ratio 
Qd/Qu = — 1/2 and of colour suppression. The ratio 
R{p^ / K*^) acquires, therefore, a tiny dependence on 
the CKM angle a. Within the SM, the numerical value 
of these NLO corrections are Ai?^ = 0.055 ± 0.13 and 
AR^ = 0.015 ± 0.11. Explicit expressions for these 
quantities, which are valid in the presence of beyond- 
the-SM physics, can be found in Refs. [54,59]. 

Further important observables are the isospin breaking 
ratio given by 



A(p7) = 



T{B+ ^ p+7) - T{B- ^ p-7) 
2(r(B0^p0^)+r(B0^p07)) 



- 1(19) 



and the CP asymmetries in the charged and neutral 
modes, 



^Cp(P7) 
Alp{pi) 



T{B- ^ p-^)-V{B+ ^ p+^) 



T{B- 
T{B° 



>p-j) + T{B+ 

P°7) - r(B° - 



P+7) 
P°7) 



r(BO ^ pO^^ + Y{BO ^ pOj) 



(20) 



(21) 



Recently, the BABAR collaboration has reported a 
significant improvement on the upper limits of the 
branching ratios for the decays B'^{B'^) ^"7 and 
p^'y. Averaged over the charge-conjugated 
modes, the current 90% C.L. upper limits are [60]: 



B{B° 

B{B^ ^ p 
B{B° iO'y) 



p°7) < 1.4 X 10" 



^7) 



< 
< 



2.3 X 10"^ 
1.2 X 10"*^ 



(22) 

(23) 
(24) 
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They have been combined, using isospin weights for 
B pj decays and assuming B{B° — > luj) — B{B^ 
^"^7), to yield the improved upper limit 



B{B fry) < 1.9 x 10" 



(25) 



Note that the equality between the p and to branching 
ratios receives S'?7(3)-breaking corrections that can be 
as large as 20%. The current measurements of the 
branching ratios for B K*^ decays by BABAR [61], 

B{B° K*"-/) = (4.23 ± 0.40 ± 0.22) x 10"^ (26) 
B{B+ -> K*+j) = (3.83 ± 0.62 ± 0.22) x 10"^ (27) 

are then used to set a 90% C.L. upper limit on the 
ratio of the branching ratios [60] : 



B[B ^ p7) 

B[B -> K*-f) 



< 0.047 



(28) 



This bound is typically a factor 2 away from the SM 
estimates [54], which we quantify more precisely in 
this letter. In beyond-the-SM scenarios, this bound 
provides a highly significant constraint on the relative 
strengths of the b d'j and b s"f transitions. 

Let us present an updated analysis of the constraints 
in the (p, fj) plane from the unitarity of the CKM ma- 
trix, including the measurements of the CP asymme- 
try a^Ks in the decays B^ / B^ J/tpRg (and re- 
lated modes) , and show the impact of the upper limit 
R{p-i/K*-i) < 0.047 [60]. The SM expressions for ex 
(CP-violating parameter in K decays), AMb^ 
B^ mass difference), AMb^ {B^-B^ mass difference) 
and a^Ks are fairly standard and can be found, for 
instance, in Ref. [62] ; the values of the theoretical pa- 
rameters and experimental measurements that we use 
are taken from Ref. [59] . The SM fit of the unitarity 
triangle is presented in Fig. 1. Note that where the 

hadronic parameters fs^ Bb^ and Cs are concerned, 
we adopt very recent lattice estimates that take into 
account uncertainties induced by the so-called chiral 
logarithms [63]. These errors are extremely asymmet- 
ric and, once taken into account, reduce sizeably the 
impact of the AMbJAMb^ lower bound on the UT 
analysis. In Fig. 1 we explicitly show what happens 
to the allowed regions once these errors are taken into 
account. The 95% C.L. contour is drawn, taking into 
account chiral logarithms uncertainties. 

As the bound from the current upper limit on 
R{p^ / K*^) is not yet competitive to the ones from ei- 
ther the measurement of AMb^ or the current bound 
on AMb^ , we use the allowed p — fj region to work out 
the SM predictions for the observables in the radiative 
B-decays described above. Taking into account these 




Figure 1. Unitary triangle fit in the SM and the result- 
ing 95% C.L. contour in the p - t) plane. The solid lines 
show the upper bounds (with and without chiral logs) due to 
AMb^, the dot-dashed lines the ones due to AMb^- The 
impact of the R{p'y / K* y) < 0.047 constraint is also shown. 
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Figure 2. The impact of the new lattice QCD estimate of 
the ratio ^ on the R{p"f / K* "f) constraint. 



errors and the uncertainties on the theoretical param- 
eters, we find the following SM expectations for the 
radiative decays: 



R^{pj/K*j) 
R°ip-f/K*-f) 



= 0.023 ±0.012, 

= 0.011 ± 0.006, 

- -0 10+°°2 

— U.lU_o.o3 J 

= -0.06 ±0.02. 



(29) 
(30) 
(31) 
(32) 
(33) 



In the CP asymmetries the uncertainties due to form- 
factors cancel out to a large extent, however, the scale 
dependence is rather large because the CP asymme- 
tries arise at the 0{as). It is interesting to work out 
the extremal values of R{pj/ K*j) compatible with 
the SM UT analysis. Any measurement of R{pj/ K*'y) 
whose central value lies in the range (0.013, 0.037) 
would be compatible with the SM, irrespective of the 
size of the experimental error. The error induced by 
the imprecise determination of the isospin breaking 
parameter ^ currently limits the possibility of having 
a very sharp impact from R{pj / K*"f) on the UT anal- 
ysis. 
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Let us comment on the impact of a recent prelimi- 
nary lattice determination of the B — > [K* , p) form 
factor at zero recoil [19]: F^' [Q) = 0.25(5)(2) and 
C = 0.91(8). In the first place, note that the cen- 
tral value of the K* form factor is in perfect agree- 
ment with the indirect determinations obtained in 
Refs. [64,65,54]. In these papers, the authors compute 
BR{B K*j) at NLL, using the inclusive channel 
B Xg-f to extract the value of the Wilson coefficient 
Cy. Their conclusion is that in order to accommodate 
the experimental data, the form factor at zero-recoil 
has to be substantially smaller than the typical light- 
cone QCD estimate (Tf^*(0) = 0.38 ± 0.06). For in- 
stance, in Ref. [54], the fitted value of the form factor 
is T^* (0) = 0.27 ± 0.04. In second place, the higher 
central value for the ratio ( given in Ref. [19] strength- 
ens the impact of the current R{p^/ K*j) upper limit. 
In Fig. 2, the additional line has been obtained using 
this new determination of the form factors ratio and 
gives a bound comparable to the Ams^ one. In order 
to fully trust this new lattice estimate of the ratio C, 
an independent cross-check of this result is mandatory; 
moreover it is necessary to analyse the old QCD sum 
rules estimate deeper so as to understand the reasons 
of this discrepancy. 

Let us finally discuss the analysis of the exclusive 
modes in supersymmetric models and entertain two 
variants of the MSSM called in the literature MFV [66] 
and Extended-MFV [67] models. In MFV models, all 
the flavour changing sources other than the CKM ma- 
trix are neglected. In this class of models there are 
essentially no additional contributions (on top of the 
SM ones) to a^Ks and AMbJAMb^, while the im- 
pact on ck, AMb^ and AMb^ is described by a single 
parameter, /, whose value depends on the parameters 
of the supersymmetric models [62]. EMFV models 
are based on the assumption that all the superpart- 
ners are heavier than 1 TeV with the exception of 
the lightest stop; no constraints are imposed on the 
off diagonal structure of the soft breaking terms. It 
can be shown [67] that under these assumptions there 
are only two new parameters in addition to the MFV 
ones, namely: S^^^i = Mlj/iMfMg) x Vtd/\Vtd\ and 
6,^i = MlJ[MiMq) X Vtll\Vts\. Where Hs the light- 
est stop mass eigenstate and A/^ is the up-squark mass 
matrix given in a basis obtained from the SCKM one 
after the diagonalization of the 2x2 stop submatrix. 
Since we are interested in the phenomenology of 6 — > c? 
transitions, we will consider here only i^^j- With 
the inclusion of this new parameter, the description 
of the UT-related observables needs one more com- 
plex parameter, g — gn + igi [67]. A signature of 
these models is the presence of a new phase in the 
B^ — B^ mixing amplitude. Using the parametriza- 
tion = r^e^'^'iMfj^, we get = [1 + / + .g| and 



A(py) 

-0 .25 
-0 . 5 
-0.75 



L = SM central values ^ = EMFV central values 
□ = SM at 68% C.L. B = EMFV at 68% C.L. 
V □ = MFV at 68% C.L., > 




SM 







0.01 0.02 0.03 0.04 0.05 0.06 0.07 0.08 



R(PY/K*y) 

Figure 3. Correlation between R{p"f / K* 'j) and A{pj) in 
the SM and m MFV and EMFV models. The light-shaded 
regions are obtained varying p, r), the supersymmetric pa- 
rameters (for the MFV and EMFV models) and using the 
central values of all the hadronic quantities. The darker 
regions show the effect of ±la variation of the hadronic 
parameters. 

9d = 1/2 arg(l + f + g). This implies new supersym- 
metric contributions to the CP asymmetry a.^Xs ■ 

The phenomenology of the MFV and EMFV models, 
analysed by scatter plots over the supersymmetric pa- 
rameter space, shows the discrimation power of exclu- 
sive modes, if one focus on ratios of exclusive observ- 
ables and their correlation. If one also scan over p 
and fj, and require that each point satisfy the bounds 
that come from direct searches, from the B — s- Xs^ 
branching ratio, and from the UT-related observables, 
one finally finds the surviving regions presented in 
Fig. 3. It shows the correlation of the isospin break- 
ing ratio A{pj) and the ratio of the branching ratios 
R{p^/ K*j). The light-shaded regions are obtained 
using the central values of the input parameters while 
the dark-shaded ones result from the inclusion of their 
la errors. In the MFV case, there are two distinct 
regions that correspond to the negative (SM-like) and 
positive C| case. For C| < 0, the allowed regions in 
MFV almost coincide with the SM ones and we do 
not draw them. For C| > 0, the allowed regions are 
different and, in general, a change of sign of both the 
CP-asymmetries (compared to the SM) is expected. 
We note that the latter scenario needs very large SUSY 
contributions to Cy , arising from the chargino-stop di- 
agrams, and for fixed values of tan/Jg it is possible to 
set an upper limit on the mass of the lightest stop 
squark. 

4 b ^ si^i~ Transitions 

4.1 B Xsi+e- 

In comparison to the B Xs^, the inclusive B — )■ 
Xgl'^i^ decay presents a complementary and also 
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more complex test of the SM. This decay is dominated 

by pcrturbativc contributions if the cc resonances that 
show up as large peaks in the dilepton invariant mass 
spectrum are removed by appropriate kinematic cuts. 
In the 'perturbative windows', namely in the low-s 
region 0.05 < s = q^/rn^ < 0.25 and also in the 
high-s region with 0.65 < s, theoretical predictions 
for the invariant mass spectrum are dominated by the 
purely perturbative contributions, and a theoretical 
precision comparable with the one reached in the de- 
cay B Xs"f is in principle possible. Regarding the 
choice of precise cuts in the dilepton mass spectrum, 
it is important that one directly compares theory and 
experiment using the same energy cuts and avoids any 
kind of extrapolation. 

In the high-s region, one should encounter the break- 
down of the heavy mass expansion at the endpoint. 
Integrated quantities are still defined; nevertheless one 
finds sizeable ^qcd/^I non-perturbative corrections 
within this region. 

The decay B XJ.^C.^ is particularly attractive be- 
cause of kinematic observables such as the invariant 
dilepton mass spectrum and the forward-backward 
(FB) asymmetry. They are usually normalized by the 
semi-leptonic decay rate in order to reduce the uncer- 
tainties due to bottom quark mass and CKM angles 
and are defined as follows: 



R\uark{s) = ^r(6 ^ X,tr)lV{h ^ X^ev), (34) 



Afb{s) 

X ^ dcosOg 



d^rjb xj+e-) 

ds d cos 6 1 



sgn(cos^£), (35) 



here 6i is the angle between and B momenta in the 
dilepton centre-of-mass frame. The so-called 'normal- 
ized' FB asymmetry, which is also often used, is given 

by 



Afb{s) 



dcosOe — — sgn(cos6'^) 

.1 ds dcosOi 

\ „ d^r{B ^ Xs£+e-) 

d cos (fe 



ds dcos( 



(36) 



For the low-s region the present partonic NNLL pre- 
diction is given by (see [68,69,72]): 



/.0.25 

io.05 



dsRluarkis) = (1.27±0.08,eaJe) X 10-^ (37) 



The error quoted in (37) reflects only the renormal- 

ization scale uncertainty and is purely perturbative. 
There is no additional problem due to the charm mass 
renormalization scheme ambiguity within the decay 
B Xsi^l~ because the charm dependence starts 
already at one loop, in contrast to the case of the de- 
cay B Xgj. The charm dependence itself leads 
to an additional imccrtainty of ~ 7.6% within the 
partonic quantity (37), if the pole mass is varied. 



pole 



/m; 



pole 
b 



0.29 ±0.02. 



The impact of the NNLL contributions is significant. 
The large matching scale nw uncertainty of 16% of the 
NLL result was removed; the low-scale uncertainty /ih 
of 13% was cut in half; and also the central value of 
the integrated low dilepton spectrum (37) was signif- 
icantly changed by more than 10% due to NNLL cor- 
rections. Using the measured scmi-lcptonic branch- 
ing ratio 6|^.p , the prediction for the corresponding 
branching ratio is given by 

B(B-.X,£+r)cut:ie[o.o5,o.25]= (38) 

^0.25 

= Selp. / ds[Rl2rkis) + Rmlis) + Rml{s)] 
J0.05 

= (1.36 ± 0.08,ca/e ) X 10-<5 

Rm^ (s) and Rm^ (s) are the non-perturbative contri- 
butions scaling with 1 / ml and 1 / ml . The recent first 
measurement of BELLE, with a rather large uncer- 
tainty [22] , is compatible with this SM prediction. 

The phenomenological impact of the NNLL contribu- 
tions on the FB asymmetry is also significant [70,71]. 
The position of the zero of the FB asymmetry, de- 
fined by ApBisa) = 0, is particularly interesting to 
determine relative sign and magnitude of the Wilson 
coefficients Cy and Cg and it is therefore extremely 
sensitive to possible new physics effects. The previous 
NLL result, where the error is determined by the scale 
dependence, is now modified by the NNLL contribu- 
tions [70,71]: 



'NLL 



0.14 ±0.02, s'^ 



NNLL 



0.162 ±0.008. (39) 



In the NNLL case the variation of the result induced 
by the scale dependence is accidentally very small 
(about ±1%) and cannot be regarded as a good es- 
timate of missing higher-order effects. Taking into 
account the separate scale variation of both Wilson 
coefficients Cg and Cy, and the charm-mass depen- 
dence, one estimates a conservative overall error on sq 
of about 5% [70]. In this s region the non-perturbative 
1 /m^ and 1 /m^ corrections to Apg are very small and 
also under control. An illustration of the shift of the 
central value and the reduced scale dependence be- 
tween NNL and NNLL expressions oi ApBis), in the 
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0.5 



0.2 0.25 



Figure 4. Comparison between NNLL and NLL results for 
Afb{s) tn the low s region. The three thick lines are the 
NNLL predictions for fi = 5 GeV (full), and ^ = 2.5 and 
10 GeV (dashed); the dotted curves are the corresponding 
NLL results. All curves for rric/mb = 0.29. 



low-s region, is presented in fig. 4. The complete effect 
of NNLL contributions to the FB asymmetry adds up 
to a 16% shift compared with the NLL result, with 
a residual error due to higher-order terms reduced at 
the 5% level. Thus, the zero of the FB asymmetry in 
the inclusive mode turns out to be one of the most 
sensitive tests for new physics beyond the SM. 

The B factories will soon provide statistics and reso- 
lution needed for the measurements of i? — > Xst^t~ 
kinematic distributions. Correspondingly, the recently 
calculated new (NNLL) contributions [68-72]'* have 
significantly improved the sensitivity of the inclusive 
B Xs£~^£~ decay in testing extensions of the SM 
in the sector of flavour dynamics. However, with the 
present experimental knowledge the decay B Xs^ 
still leads to the most restrictive constraints as was 
found in [74]. Especially, the MFV scenarios are al- 
ready highly constrained and only small deviations to 
the SM rates and distributions are possible; therefore 
no useful additional bounds from the semi-leptonic 
modes beyond what are already known from the B 
Xs^ can be deduced for the MFV models at the mo- 
ment. Within the model-independent analysis, the im- 
pact of the partial NNLL contributions on the allowed 
ranges for the Wilson coefficients was already found 
to be significant. In this analysis, however, only the 
integrated branching ratios were used to derive con- 
straints. It is clear that one needs measurements of 
the kinematic distributions of the B Xsl~^£~, the 
dilepton mass spectrum and the FB asymmetry in or- 
der to determine the exact values and signs of the Wil- 



-0.5 




*We add here that the three-loop mixing is fully are under con- 
trol. A quite recent calculation of a missing NNLL mixing piece 
leads to a correction below 2% [73] 



Figure 5. Four different shapes of the normalized FB 
asymmetry Afb for the decay B — > Xs£'^£~. The four 
curves correspond to four sample points of the Wilson coef- 
ficients that are compatible with the present measurements 
of the integrated branching ratios. 



son coefficients. In fig. 5, the impact of these future 
measurements is illustrated. It shows the shape of the 
FB asymmetry for the SM and three additional sample 
points, which are all still allowed by the present mea- 
surements of the branching ratios; thus, even rather 
rough measurements of the FB asymmetry will either 
rule out large parts of the parameter space of extended 
models or show clear evidence for new physics beyond 
the SM. 



4.2 B -> K^*k+£- 

For completeness, let us briefly comment on the im- 
pact of exclusive B K^*^£^£~ modes. First of all, 
let us stress that hadronic uncertainties on these exclu- 
sive rates are dominated by the errors on form factors 
and are much larger than in the corresponding inclu- 
sive decays. In fact, following the analysis presented 
in Ref. [74], we see that inclusive modes already put 
much stronger constraints on the various Wilson coef- 
ficients. 

Concerning the measurement of a zero in the spec- 
trum of the forward-backward asymmetry, things are 
different. According to Refs. [75,64] the value of the 
dilepton invariant mass (qq), for which the differential 
forward-backward asymmetry vanishes, can be pre- 
dicted in quite a clean way. In the QCD factorization 
approach at leading order in Aqcd/^6, the value of 
Qq is free from hadronic uncertainties at order (a 
dependence on the soft form factor and the light 
cone wave functions of the B and K* mesons appear 
at NLL). Within the SM, the authors of Ref. [64] find: 
Qq = (4.1 ± 0.6) GeV^. As in the inclusive case, such a 
measurement will have a huge phenomenological im- 
pact. 
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